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A theory of coherent multiple Andreev reflections (MAR) is developed for superconductor-normal 
metal-superconductor (SNS) interferometers. We consider a Y-shape normal electron beam splitter 
connecting two superconducting reservoirs, where the two connection points to the same supercon- 
ductor can have different phase. The current is calculated in the quantum transport regime as a 
function of applied voltage and phase difference, I{y,(f)). MAR in interferometers incorporates two 
features: interference in the arms of the splitter, and interplay with Andreev resonances. The latter 
feature yields enhancement of the subgap current and current peaks with phase-dependent positions 
and magnitudes. The interference effect leads to suppression of the subgap current and complete 
disappearance of the current peaks at = n. The excess current at large voltage decreases and 
changes sign with increasing phase difference. 
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I. INTRODUCTION 

Mesoscopic circuits with multiterminal electron wave 
guides have interesting and useful physical properties. 
Figure |l| shows an example of a circuit where a Y-branch 
wave guide acts as a coherent beam splitter for electrons 
injected from a normal or superconducting reservoir. 
Such a circuit with normal electron reservoirs (NYN) 
has been first suggested to test statistical properties of 
the Fermi electrons by measuring |Current-current corre- 
lations in the arms of the splitter E'0 This correlation is 
negative for the case of FermLstatistics, which was indeed 
observed in the experiment .0 However, the correlations 
may change the sign and become positive (similar to the 
case of Bose particles) when the electroas are injected 
from a superconducting reservoir (NYS).Q This effect is 
related to the Cooper pairing in the superconductor. 

Another kind of interference effect, narr^f a phase de- 
pendent conductance, has been predictedBQ for a circuit 
geometry (SYN) where the arms of the electron beam 
splitter are connected to a superconductor. In this case, 
the electrons undergo Andreev reflection from the NS 
interfaceaJ, picking up the superconducting phase at the 
connection points. The phase difference at the connec- 
tion points is created by a supercurrent flowing along 
the surface of superconducting electrode, e.g. due to a 
presence of a magnetic field (see Fig. ^. The phase de- 
pendence of the conductance of such NS interferomet^s 
has been observed in a large number of experiments.ETtEl 




FIG. 1. Schematic picture of the SNS interferometer, which 
consists of the superconducting loop, to the left, with two 
electrodes and 52 that are connected via the Y-branch 
wave guide to the right injection electrode 5*3, which can be 
superconducting or normal. The two left electrodes in the 
loop and the Y-branch wave guide constitute an SNS junc- 
tion (S1NS2), with bound Andreev states broadened due to 
the coupling to the injection electrode S3. The phase differ- 
ence over the S1NS2 junction is controlled by means of the 
magnetic flux through the loop. 

In this paper, we investigate the properties of current- 
voltage characteristics of SNS interferometers, i.e. SYS- 



circuits where the arms of the Y-branch splitter and 
the injection lead are connected to superconducting elec- 
trodes, as shown in Fig. |l|. As is well known, the cur- 
rent between two superconductors at applied voltage 
eV < 2A (A is the superconducting gap) is governed h« 
the mechanism of multiple Andreev reflections (MAR).U 
The fingerprint of this transport mechanism is t>ip snb- 
harmonic gap structure (SGS) at eV = 2A/n.tM3 In 
SNS interferometers one should expect the SGS to be 
significantly modified and become sensitive to the super- 
conducting phase difference (see Fig. |l|). This conclu- 
sion is supported by the following qualitative argument. 
One may consider the interferometer as an SNS junction 
{S1NS2) which is coupled to the injection electrode S3 
via the Y-branch wave guide. The SNS junction con- 
tains bound Andreev levels which are broadened due to 
the open connection to the injection lead. The Andreev 
levels, the energy of which depends on the superconduct- 
ing phase difference, produce resonances in the current 
transport between the superconducting reservoirs. These 
resonances will significantly modify the SGS, similar to 
the effect |-0|f superconducting bound states in long SNS 
junctionsJlil However, in the SNS interferonoeters the res- 
onance structures will be phase dependent Jlj 

The problem has become particularly interesting due 
to recent experiments on a thxee-terminal SNS interfer- 
ometer by Kutchinsky et al.,t3 in which phase depen- 
dence of the conductance structures has been observed. 
In this experiment, the normal region of the interferome- 
ter was fabricated with diffusive two-dimensional metallic 
film. The effect will be more pronounced with a ballis- 
tic normal region, especially in the quantum transport 
regime, when the Y-wave guide supports a small amount 
of conducting modes. In practice, such devices could be 
fabricated with etched or gated ballistic 2D [electron gas 
in a multilayered semiconducting structure.tZl Quantum 
transport in NS interferometers has been investigated in 
Refs. p^|-pO|, and the resonance enhancement of the two- 
particle Andreev current has been discussed. In SNS in- 
terferometers the situation is more complex, since the in- 
terplay between the Andreev bound states and MAR will 
lead to resonances appearing in all multiparticle currents. 
The purpose of our study is to develop MAR theory for 
SNS interferometers, i.e., to include the interference ef- 
fect and the effectj-of Andreev resonances in the MAR 
calculation schemeE3 and to analyze the phase dependent 
current- voltage characteristics (CVC). 

The paper is structured as follows. In section II we 
present the model of the interferometer and derive equa- 
tions for the MAR scattering amplitudes. In section III, 
we present the results of numerical calculation of the 
CVC of the interferometer, and in section IV, we present 
analytical perturbative analysis of the resonant current 
structures in the CVC. 
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II. SCATTERING STATES 

We consider the SNS interferometer shown in Fig. 
which consists of a normal-electron Y-branch quantum 
wire, fabricated e.g. with etched or gated 2D electron 
gas, which is connected to two bulk superconducting elec- 
trodes. The superconducting electrode connected to the 
two arms of the splitter is ring shaped, so that magnetic 
flux can be sent through the ring to induce a supercon- 
ducting phase difference </> between the connection points 
Si and 82- The current is sent through the interferome- 
ter by applying voltage V between the electrode 6*3 and 
the ring. 

Our aim is to calculate the injected current as a func- 
tion of the applied voltage and the phase difference, 
I{V, 4>). To this end we shall apply the theory of coher- 
ent multiple Andreev reflections (MAR), which is based 
on the calculation of scattering states created by in- 
coming quasippiticles from the three superconducting 
terminals Ej'EJa For the sake of simplicity, we adopt a 
single mode description of the wires and model the split- 
ter by a 3x3 scattering matrix S (see Eq. (^) below). 

Solving a coherent MAR problem includes the follow- 
ing steps: First the recurrence for the scattering ampli- 
tudes is formulated by employing the boundary condi- 
tions at the scatterer and at the NS interfaces. Then 
a solution for the recurrence can be constructed, from 
which the current is calculated. In the present case of 
a three-terminal junction standard MAR technique can- 
not be used directly. To solve this problem we will de- 
scribe the Andreev reflections from the two interfaces of 
the interferometer, seen by the injector S'3, as a com- 
plex reflection from a single effective interface; by doing 
this we reduce the 3-terminal problem to an equivalent 2- 
terminal problem, on which we can apply standard MAR 
technique for a 2-terminal junctions developed in Ref. |2^. 
Reflection from the interferometer (left interface) is then 
represented by a non-trivial transfer matrix that contains 
information about Andreev states in the interferometer. 
To more clearly understand the effect of the phase differ- 
ence on the injection current, we will neglect the length 
of the normal wires and assume all NS interfaces to be 
completely transparent. This will remove unnecessary 
complications due to normaLelectron resonancesEj and 
superconducting resonanceslij at^ciated with the Saint- 
James - de Gennes bound states.tj At the same time, the 
resonant property of the junction as well as the phase 
sensitivity of the current will persist due to the pres- 
ence of Andreev states in the interferometer with energy 
E = ±E„„ 



Ea - AJl- Dsm^{(b/2). 



(1) 




FIG. 2. A close up of the Y-branch wave guide in the SNS 
interferometer (upper figure) and in the effective 2-terminal 
junction where the left electrode has phase dependent scatter- 
ing properties (lower figure). We consider a symmetric phase 
drop over the SNS junction (the electrodes 5*1 and 5*2) and 
a symmetric voltage drop eV from the loop to the injection 
electrode S3. The directions of the scattering amplitudes a„ 
and /3,i (for electrons) and 7„ and S„ (for holes) are shown 
with arrows. 

In the MAR regime, the time-dependent scattering 
states consist of a superposition of scattered waves with 
energies En = E + neV {n is an integer) shifted with 
respect to the energy of the incoming wave E. In each 
branch of the wire the wave functions of the sideband 
n consists of a superposition of electron and hole plane 
waves moving in both directions, the corresponding am- 
plitudes being labeled with a„, ...(5„ according to Fig. ||. 
In the normal wires, the scattering state wave function 
has the form, 
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where ^'^(i) is the wave function in the i-th wire at the 
left side {i = 1,2), x is a coordinate along the wire, 
^'^ is the wave function at the right side and — 

^j2m{EF ± En) is the normal wave vector {% = \). In 
the superconducting electrodes, the corresponding wave 
functions, including six possible source terms, are: 
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where An, ■ 

elemental 

equatioECJ, 

where 



..Ln are constant coefficients and ip^ 
solutions to the Bogohubov-dc Gcnncs 
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CT„ = sign(£;„), fc„ = yj2m{EF ± Cr„Cn): 

and is the Pauh matrix. The different source terms 
are distinguished by the index j — {1, 2, 3, 4, 5, 6} corre- 
sponding to injection of an electron/hole from each of the 
electrodes. Injections from the left and from the right 
generate essentially different scattering states and cur- 
rents, and to keep track of the side of injection we will 
introduce an additional side index z/ = =pl corresponding 
to injection from the left electrodes and right electrode, 
respectively. 



v = ^l ^ jG {1,2,3,4} 
V = 1 <J=^ j e {5, 6} 
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The scattering in the Y-branch splitter is described by 
a unitary scattering matrix: 
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where the scattering matrix, 
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is chosen to be symmetric, 6*31 = ^32 = t, with a real 
transmission amplitude from right to left, t > 0, and 
also to be independent of energy and thus the same for 
electrons and holes. 



(10) 



The latter approximation is reasonable since the scale of 
the energy dispersion of the scattering matrix is given 
by the Fermi energy, while we are interested in a much 
smaller energy interval comparable to the superconduct- 
ing gap, Ep ^ A. The unitary conditions for S are 
expressed through the equations, 1 = R + D + T, 1 — 
|rop -1- 2T, where D = R = jrp and T = and 



-(r + d)*, r-d 



(11) 



The transmission coefficient D concerns the transparency 
of the SNS junction of the interferometer and controls 
the position of the Andreev levels, Eq. (Q), while the 
transmission coefficient T describes coupling of the SNS 
junction to the injection lead. 

By means of a canonical transformation the scattering 
between the injection electrode and the interferometer 
can be separated from the scattering between the arms 
of the interferometer. To this end we introduce new scat- 
tering amplitudes, = (ai,n ^ o;2,rt)/V2, and similarly 
for /3„,7„, and (5„, and rewrite Eq. (0) on the form, 



S' = 



S' 



2T 
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Pn ={r~ d)a;_ 



(12) 



(13) 



The scattering equations for holes have a similarly form. 
It is clear from Eq. ( |l2|) that the transport through the 
interferometer, from the right electrode to the left elec- 
trodes, is independent of the (-)-coefficients. Thus, we 
can treat the junction as an effective two-terminal junc- 
tion with transparency 2T. 

It is convenient for the further calculations to introduce 
vector notations. 



Pi 



(14) 



and similar for the other scattering amplitudes, and to 
rewrite Eq. (fl^) through a transfer matrix, 



Ta' 



(15) 



where T = (1/V2r)(l - {role-'^^Pa^). The same equa- 
tion also holds for the hole coefficients. It is convenient to 
gauge out the reflection phase p from the transfer matrix, 
which can be done by a transformation of the scattering 



amplitudes, q+ 



whereby 



T^(l/V2T)(l-|ro|a.). 



(16) 



We now proceed with a derivation of the recurrence 
relation and first consider the right NS-interface, where 
we match the wave functions in Eq. (H) and (H), 



7« = K +i;,-(5i^(5on- 



(17) 



This is the recurrence relation for an ideal SN-interface 
where the Andreev reflection amplitude is given by 



and the source term. 



(18) 
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corresponds to the injection from the right electrode, 
u = 1. Similar relations hold for each of the two SN- 
interfaces at the left side. In terms of the vectors and 
^ introduced in Eq. (M), they have the form, 
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for n 0. For n — there is also a source term describing 
injection from the left electrodes = —1). Now the 
vectors and 7^ can be eliminated by means of Eq. 
( |l3| ) for electrons and for holes, which yields 



(21) 



where 



cos ■ 



(22) 



6 = 6' + p. The source term, Yj , for j = 1,2 has the form 
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where 6 = ii7o(A/2^) tan(0/2). The source term for the 
injection cases j = 3, 4 has form similar to Eq. ( p3[ ) with 

The matrix 1/^ describes complete reflection from the 
interferometer for energies inside the energy gap, \En\ < 
A, where it obeys a standard transfer matrix equation. 



K^fT.l/t^ = a,. For 



this is purely Andreev reflec- 



tion, identical to the one described by Eq. (|l7|), which 
implies that in this case the interferometer works as a 
single ideal SN interface. In the general case, 4> ^ 0, the 
reflection consists of both Andreev and normal reflec- 
tions, the probability of Andreev reflection being given 
by the matrix element |(T^„^)ii|~^ ~ cos^((/)/2). Thus, 
at = TT the probability turns to zero, and therefore the 
Andreev transport through the interferometer is blocked. 
Furthermore, the transfer matrix V^^ contains informa- 
tion about Andreev bound states in the interferometer. 
Assuming for a moment the SNS junction disconnected 
from the injection lead, T = and considering a sta- 



tionary version of Eq. (|T]), at = 7^, in combination 



with Eq. (|T^) we get the solvability condition for these 
equations on the form Im[(y^)ii — (1^)12] = 0, which 
gives the Andreev level in Eq . (p| ). 

In principle the equations (|17|), ( 2l) ) and (15) provide 
a complete set of recurrences for the MAR amplitudes. 



However, following Ref. |2^, it is convenient to introduce 
new amplitudes, c„±, which allow us to get rid of redun- 
dant MAR amplitudes and to unify notation for different 
injection cases, 



C2m+ = + 72m'^l'' 
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Then the recurrence relations Eq. (|17|), ( |2l| ) and (|15D can 
be written on a compact form. 
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(25) 
(26) 

(27) 
(28) 



and the index fi — {—l)"^ specifies the particular form of 
the transfer matrices Un and T„ for different side bands 
and for different injection directions. 

Equations (p4[)-('28|) realize a mapping of the MAR 
problem on the problem of wave propagation along the 
energy axis. The wave amplitudes for the propagation 
in the upward and downward directions are given by the 

T I 

upper, c„_|_, and lower, c„_|_, components of the vector 



(29) 



The probability current flowing along the energy axis is 
defined in the usual way as 



Jn± 



-nil 



(30) 



This probability current is conserved within the super- 
conducting gap, j^j. =const for \En\ < A, due to the 
unitarity properties of the matrices r„ and [/„ . Violation 
of the unitarity for Un outside the gap indicates leakage 
of the probability current into the reservoirs. We notice 
that in real space, the probability current is carried al- 
ternatively by electrons and holes. 

The charge current through the interferometer is de- 
termined by the probability currents of all the side bands, 



ME) 



2tt 
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dE J^{E)nF{E), 



(31) 



(32) 



j/,n— even 



In this equation, the spectral current Jv{E) consists of 
the sum over the scattering states and all the sidebands 
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at the injection side of the junction, v\ the factor \E\/£_ is 
the superconducting density of states, and np{E) is the 




FIG. 3. Schematic picture of the MAR ladder in energy 
space, for an injected quasiparticle at energy E. The inter- 
faces of the junction are shown as vertical lines, where at 
the right side there is an ordinary NS interface, while at the 
left side (double line) there is the effective phase dependent 
SN-interface. The effective scatterer is shown as a dashed line 
in the center of the junction. The two-component scattering 
amplitudes c„± represent up and down going electrons (solid 
line) and holes (dashed line). The corresponding probability 
current at level n is j^j^. 

For n > 0, the free solution of Eqs. ( p5| ) and ( p6| ) can 
be written on the form, 

M„o = fn-lUn-lfn-2-Ulfo, (33) 

while the solution for negative n, n < 0, reads 

c„+ = M~o^co_ ,c„_ = C/,7^c„+. 

At the injection point n = 0, the MAR coefficients are 
found from Eq. (25), 



UoCa-+Yi. 



(34) 



This equation imposes two relations between the four un- 

T I 

known coefficients, Cg^, Cg^; another two relations are 
imposed by the boundary conditions at infinity in energy 
space 



hnr c„± = 0, 

n— >±oo 



(35) 



which fix the ratios Cg^/cQ_|_. To implement the boundary 
conditions, we introduce effective reflection coefficients, 
r„± , which characterize the intensity of the back scatter- 
ing at specific points of the MAR ladder. The definition 
reads, 



= c,Vr„+, n > 
= cl^iayf^-, n < 0, 



where 
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(36) 



(37) 



With help of the vectors f„± the solution in Eq. ( [3^ ) 
can be presented in a compact form (using the equality 
det(M„™) = 1), 



c„+ = r„ . , n > 0, 



ZnO 



(38) 



where the brackets denote a scalar product of two- 
component vectors, and 

ZnO = {f*o^,U^'MjU-'r^+). (39) 
The corresponding expression for negative n is 



(Tyfn-, n < 0, 



(40) 



ZOn = {ro+,<TyUoMonUnCryrn-)- 



The form of solution in Eqs. (ps])- (^) is particularly 
useful for the calculation of the SGS: The essential pro- 
cesses which contribute to the SGS involve transmission 
through the energy gap region. This transmission con- 
tains resonances which are included in the transfer matrix 
M„g. The latter is straightforward to calculate because 
it does not require excursions to infinity. On the other 
hand, the calculation of the quantities r„± includes such 
excursions which, however, go outside the energy gap, 
and the corresponding recurrences rapidly converge. 

To calculate the SGS, it is convenient to separate out 
the current associated with the scattering across the gap, 
which is most important for the SGS, from the current 
of thermal excitations which involves transitions between 
states below or above the gap. To this end, we first 



rewrite the spectral current in Eq. (32) through leakage 



currents into the electrodes defined as, j„ 



By using the equality 
d^), we get 



■P 

Jn-- 



following from Eq. 
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(41) 



(here we explicitly write the scattering state index j, 
introduced in Eq. (||)). We note that the conserva- 
tion of the probability current inside the energy gap 
implies that the leakage current is zero inside the gap, 
jn,] = 0, \En\ < A. Having made this observation, we 
write the total current in Eq. (p2[) on the form. 
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(42) 



The first two terms in this equation correspond to the 
current of thermal excitations: they cancel each other 
when the, temperature approaches zero due to detailed 
balance, E3 



for the partial currents defined as 



(43) 



(44) 



Keeping the last term in Eq. (|42|), we get for the current 
at low temperature, T <C A, the following equation, 



/ = ^%yn-2A)(/„^+/^), 
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(45) 
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The current Jn,v{E) in Eq. ( p4[ ) can be presented after 
some algebra on the form which is useful for analytical 
study, 



{E)inAE) 
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where 



in.i^{E) = (f„+, [f/t^V^J/,^^ - cr^]r„+). 



(47) 



(48) 



III. NUMERICAL IVC 

In this section, we present the results of numerical cal- 
culation of the current given by Eq. (^ ) , where the spec- 
tral current is given by Eq. ( p7| ) . In order to numerically 
evaluate the effective reflection amplitudes r„± we use 
the equations, the definitions Eq. (^), using Eq. 
and Eq. (IsF 



rn+ 



lim 



M, 



t(22) 
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lim 



M, 
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M, 



(11) 



(49) 



which follow from the definitions in Eq. (|36|), and also 
from Eq. (^) and Eq. (|35|). For zero phase difference, 
the current-voltage characteristic (CVC) in Fig. || shows 
the well known subharmonic gap structure for short junc- 
tions with current structures at eV = 2A/n, cmpei- 
sponding to thresholds for the n-particle currentsEIlEJ. 
These current structups are associated with the usual 
gap edge singularities ,c3 and they are most pronounced 
if the transparency 2T of the junction is small, T ^ 1. 



The CVCs of the interferometer for different values of 
the phase difference at zero temperature are shown in 




1 1.5 2 2.5 

eV/A 

FIG. 4. The subharmonic gap structure of a supercon- 
ducting point contact with current steps a eV = 2A/n, 
n = 1, 2, 3, ... (most visible on the logarithmic diagram). The 
transmission probability of the point contact is Tpc = 0.08, 
which corresponds to an SNS interferometer with (p — and 
2T = 0.08. 

Increasing the phase difference, the structure changes 
and the positions of the current structures can no longer 
be associated with the subharmonics of the energy gap. 
Instead, we observe enhancement of the subgap current 
and the appearance of large current peaks which move 
downwards in voltage with increasing phase difference. If 
the transparency is small, as in Fig. 0, the peak positions 
reflect the Andreev spectrum of Eq. ffl) : the peaks appear 
at eV ~ Ea and eV w A + Ea, as shown in Fig. ||. 
With further increase of the phase difference, the current 
structures decrease and completely disappear at (j) = n. 
At the same time, the excess current at large voltage 
becomes large and negative, turning into a deficit current. 
The CVC is 27r-Deriodic and svmmctric around tt. 




0.4 

(t)/jl 0.8 



eV/A 
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FIG. 5. The current voltage characteristics at zero temper- 
ature of an interferometer with 2T = 0.08, 7? = 0.1 and (j) = Q, 
(j) = 27r/5, (f) = 3n/5, 4> = O.Ott and = 0.997r. 




FIG. 6. Left: the energy of the (positive) Andreev level for 
D = 0.9, as a function of the phase difference. Right: the 
current peak (light region) at eV « A + Eaifp) follows the 
position of the Andreev level. 



It is illuminating to plot the individual n-particle cur- 
rents, given by Eq. (^6|). In Fig. ^ the first four n-particle 
currents are shown for (f) — St:/5 and for the same inter- 
ferometer parameters as in Fig. |[ Higher order currents 
are too small to be visible in the same picture, and their 
contribution to the current may be neglected. For com- 
parison, the inset shows the corresponding CVC when 
(j) = 0. It follows from the plots for (f> ^ 0, that the 
peak at low voltage originates exclusively from the four- 
particle current. It exists within the voltage interval 
A/2 < eV < A and within a certain window of phase 
difference. The next peak, at larger voltage, results from 
a combination of two structures: the peak of the three- 
particle current, and the overshoot of the onset of the 
pair current. The onset of a single-particle current at 
eV > 2A is reduced in comparison with the ip — case, 
which leads to significant reduction of the current at large 
voltage. The reason is spectral intensity transfer from the 
continuum to bound Andreev states. The large deficit 
current at high voltage is apparently the result of this 
reduction as well as of the above mentioned reduction of 
the multiparticle currents related to the suppression of 
Andreev reflections by the interferometer. The further 
analysis shows that the single- and three-particle cur- 
rents (and in general odd n-particle currents) injected 
from the electrode 5*3 and from the electrodes 5*1^2 are 
identical. In contrast, the two- and four-particle currents 
possess large structures only for injection from the elec- 
trode 53 while currents injected from the interferometer 
electrodes 51^2 are negligible small. In the next section, 
we will proceed with the interpretation of these features. 




eV/A 

FIG. 7. The current voltage characteristics of the same in- 
terferometer as in Fig. |^ with (j) = 3n/5. The total current is 
shown with dotted line. The single-particle current (a) is sup- 
pressed compared to the case when = 0. The two-particle 
current (b) of quasiparticles injected from the right electrode 
has a significant onset at eV — A + Ea, while the two-particle 
current (c) injected from the left is small. The three-particle 
current (d) has a peak at eV ~ A + Ea] the four-particle 
current (e) has a peak at eV « Ea- The inset shows the 
corresponding n-particle CVCs for = 0. 

IV. RESONANCE APPROXIMATION 

In the case of weak coupling of the S1NS2 junction 
in the interferometer to the injection lead S3, T <^ 1, 
it is possible to perform a perturbative analysis of the 
n-particle currents in order to explain the origin of the 
phase dependent current structures in the CVC. In this 
limit, the current is dominated by the contribution of the 
Andreev resonances, which we will consider in the lowest 
order with respect to the small parameter T ^ 1. We 
will proceed with a detailed analysis of the four lowest- 
order n-particle currents in order to explain the major 
structures in the IVC. 

Let us start the discussion with the approximation for 
the effective refiection amplitudes at specific levels in the 
MAR ladder, f„±, in Eq. (^7|). In the lowest order ap- 
proximation with respect to T, these reflection ampli- 
tudes are contributed by single back scattering. 

Since the denominator |2;„oP in expression Eq. (^) for 
the spectral current contains a factor l/(2r)" which 
comes from the product of transfer matrices Tn in Eq. 
(^), we can neglect all terms except the zero order one 
in Eq. (^, whereby 

_ 4\E^\{El-El) 4|g„|g„ 

A2 ^„cos2| ' A2 ■ ^"'^ 
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The denominator |z„oP in Eq. has to be consid- 
ered separately for different currents, n, and injection 
cases, V. 

The single-particle current involves only a single cross- 
ing of the tunnel barrier, schematically shown in Fig. 

Thus, it cannot be resonant, and the denominator 
|ziop can be considered in the lowest-order approxima- 
tion. The expression for the single particle current, de- 
rived from Eqs. (Hq), (|47|), and (M) then reads. 



h 



2Te 



-A 



dE NLiE)NR{E + eV), 



A-eV 



where 



NUE) 
Nr{E) 



\E\VW^^ 

- El 
\E\ 



VE^ 



(52) 

(53) 
(54) 



The form of Eq. (^2|) is a standard one for a tunnel cur- 
rent given by the tunnel Hamiltonian model were Nf({E) 
is the density of states (DOS) in the injection lead, while 
iVi (E) has the meaning of the effective phase dependent 
DOS in the interferometer. An important property of the 
DOS of the interferometer is the disappearance of the sin- 
gularities at the gap edges for nonzero phase difference. 
The reason is that the resonances (Andreev states) move 
down into the gap, according to Eq. (|l|) (see Fig. ||). 
Consequently, the onset of the single-particle current at 
eV — 2A, which is of the order T when = 0, is reduced 
at finite phase difference (see Fig. |^). 



-1 



-2 



-3 





2 10 12 

Nl Nr 

FIG. 8. The effective superconducting density of states in 
the left and the right electrodes of the SNS interferometer, 
with D = 0.9 and (ji = 3tv/5. The shaded regions represent the 
states filled at zero temperature. When the phase difference 
is changed from zero, the bound Andreev states split from the 
continuum, and the peaks in the density of states at the gap 
edges are smeared in the left electrode. 



In order to calculate the pair current we will first 
inspect the corresponding transfer matrix, M^q^ = 
Tf^^tj^^T^^ . For V — 1 (injection from the right) this 
transfer matrix has zeros in the limit T = at the ener- 
gies of the Andreev bound states in the interferometer. 
El ~ ±Ea. The expansion of the matrix in powers of T 
has the form 



1 



T|Ci|cosf 



[{Ef~El)ia.^~^ay){l-T) 



+ iT{Ei\^i\ - a^VDRsin^ |)] + 0{T^), (55) 

where we have used the identity, sin 9 ~ derived 
from the unitarity condition Eq. ( [ll| ) and 9 = 9' -\- p, m 
the limit T — Q. Equation ( |55| ) clearly demonstrates the 
resonant behavior. The first term is real and proportional 
to the deviation of the energy Ei for the Andreev reflec- 
tion from the energy of the Andreev state. This term 
determines the position of the resonance. The second, 
imaginary, term is proportional to T and determines the 
width of the resonance. Both the position and width of 
the resonance depend on the phase difference. At T <C 1 
the spectral current is dominated by the Andreev reso- 
nance, and neglecting the non-resonant current we arrive 
at the expression 



^2 



2e 



dE^ 



r+r_ 



A-eV 



[El - Eaf + C-^^Y'' 



where 



T±=TD±NR{EaTeV), 



E-.± 



"^"^^ 2 



(56) 

(57) 
(58) 



< D± < ^/D. a factor of 2 appears in the current 
expression, in addition to the two-particle factor n = 2 
in Eq. (^7|), since there are two Andreev resonances, at 
Ea and —Ea- The integrand in Eq. ( |5^ ) has the form of 
the transmission coefficient for an asymmetric quantum 
mechanical double barrier structure. The difference in 
the tunnel rates r± is due to the different DOS for the 
injection and exit energies, but also due to the different 
effective transparencies of the Y-branch splitter, TD±, 
for electrons and holes. The resonant transparency of the 
junction for the pair current is proportional, according to 
Eq. (^), to the product 



A2(A' 



^eos^^ 

4K2 cos 2, 



(59) 



which turns to zero at = tt. Thus the resonant pair 
current is blocked at = tt, in agreement with our ear- 
lier general conclusion about the suppression of Andreev 
refiection at the interferometer in this case. Evaluation of 
the integral over energy in Eq. ( p6[ ) under the condition, 
Ea/ A < 1 — T, yields the pair current 
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R_ 4er+r. 



Tit _ 
^9 — 



(60) 



The magnitude of this current is generally of order T; 
however, in the vicinity of the resonant onset, eV — 
A + Ea + eSV for deviations eSV, the rate r+ increases 
due to the enhancement of the DOS in the right super- 
conductor, r+ ^ which yields the overshoot seen 
in Fig. 1^. At lower voltage, A < eV < A + Ea, the pair 
current is non-resonant and proportional to T^, which is 
also the case for the pair current I2 injected from the 















Ea' 















FIG. 9. Schematic pictures of important scattering pro- 
cesses. Resonances are shown as circles, a; Non-resonant 
single particle current; b: resonant pair current from the right 
side involves one resonant Andreev reflection, two equivalent 
resonant paths are possible; c: resonant 3-particle current, 
involves one resonant Andreev reflection and one reflection 
by the quasi-resonance at the energy gap edge, an equivalent 
path exist also for injection from the left; d: resonant 4-parti- 
cle current involves two resonant Andreev reflections (double 
resonance) . 

The resonant three-particle current can be analyzed in 
a similar manner. We construct the corresponding trans- 
fer matrix using the pair current matrix, Eq. (|55|), and 
for ^ = 1 we get M^^^ = M2q^U2^T2^. This equation 



has resonances at the same energies as the pair current, 
El = ±Ea, but nevertheless the resonant current is small, 
/g ~ T^, since the MAR path, Fig. includes an ad- 
ditional crossing of the tunnel barrier. The magnitude 
of the current is considerably enhanced in the vicinity of 
voltage eV = Ea + A when the second Andreev reflection 
occurs at the gap edge of the right electrode, E2 = A (see 
Fig. ||c). The singular DOS at the gap edge in this elec- 
trode appears in the MAR calculation as a quasibound 
state situated at the gap edge. Thus, this MAR process 
can be interpreted as an overlap of two resonances, one at 
the Andreev bound state at energy —Ea, and another at 
the quasibound state at the gap edge. A similar analysis 
applies to the current injected from the interferometer. 



ly = —1: in this case, the double resonance is formed 
by the Andreev state at Ea and the quasibound state at 
lower gap edge. The height of the double resonance peak 

A double resonance involving two Andreev bound 
states occurs in the four-particle current injected from 
the right as shown in Fig. ||d. The corresponding MAR 
path is selected by the conditions Ei = —Ea and E^ = 
Ea- The resulting peak is situated at eV = Ea and has 
the hight {l4)max ~ T. In the vicinity of the peak, 
e6V — eV — i?a <C A, the expression for the resonant 
four-particle current becomes. 



r« — __ 

^4 " 27r 



dEr 



where 



{e6VY-\{Vl+Tl) 



TD^NR{2Ea). 



Ein 



(61) 



(62) 



The current peak under consideration is a full-scale cur- 
rent structure in a region of the IVC which in point con- 
tacts is dominated by non-resonant four-particle current, 
which is of the order T^. As any Andreev resonance, the 
four-particle current peak strongly depends on the phase 
difference; it disappears at = tt. The peak exists within 
the voltage interval A/2 < eV < A and within the phase 
difference window < sin^ (j)/2 < 3/ (4D). 

We note that the positions of all discussed current 
structures are proportional to the order parameter A and 
therefore temperature dependent and scale with A. 

We end this section with the discussion of CVC at large 
voltage, eV 2A. The current at large voltage results 
exclusively from the single particle and pair currents. It 
is straightforward to calculate the single particle current 
at large voltage from Eq. (52), 



2T 



eV-A 



dEiEi^iEp _ 2Te^V 
{El-El)^,~^^ 



(63) 



The first term in this equation is the current through 
the normal junction, while the second term is the single- 
particle contribution to the excess current. 



h.. 



eTVDA\sm{^/2)\ + 0{T^ 



(64) 



This excess current is large, ~ T, and negative, which 
reflects the suppression of the current onset at eV = 2A. 
The pair current contribution to the excess current is 
dominated by the resonance and is given by Eq. (|6^). 
The total excess current to leading order in T has the 
form 



E. 



-eTAVDR 



|sin(0/2)|= 



1 - £)sin^(</>/2) 



(65) 



The excess current in the interferometer is large and neg- 
ative, in sharp contrast to quantum point contacts, where 
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the excess current is small, ~ T^, and positive .c3 How- 
ever, it turns to zero at = 0, in agreement with our 
observation that in this case the interferometer behaves 
as a quantum point contact. The excess current achieves 
its minimum value at 6 = tt because of blockade of the 



< 




FIG. 10. The first order (~ T) excess current as a function 
of the phase difference for different transparencies D of the 
SNS junction of the interferometer. 



V. CONCLUSIONS 

We have investigated the phase dependence of the 
dc current-voltage characteristics (CVC) in SNS inter- 
ferometers consisting of two superconducting reservoirs 
connected by a Y-shaped normally conducting quantum 
wire. We have developed a generalization of the exist- 
ing coherent MAR theory, which fully incorporates the 
effect of interference, as well as resonance effects due to 
Andreev bound states in the interferometer. 

The phase dependent current-voltage characteristic 
I(V, (f) has been numerically investigated for single-mode 
wires with length smaller than the superconducting co- 
herence length. The analytical study of the current was 
performed using a resonance approximation in the limit 
of weak coupling of the interferometer to the injection 
lead. 

We found significant enhancement of the subgap cur- 
rent at finite phase difference compared to the case of 
zero phase difference. The enhancement is accompanied 
by a shift of the onset of the current at 2A to smaller 
voltages, A -I- Ea[(j)) < 2A, and appearance of current 
peaks at Ea and A + Ea- These features are produced by 
resonant Andreev reflection of the injected particles by 
the S1NS2 junction of the interferometer: they are most 
pronounced for weak coupling of the S1NS2 junction to 
the injection lead. The positions of all current structures 
are temperature dependent and scale with A. 

We further found strong suppression of the current 



when the phase difference approaches (j} —f tt: in this re- 
gion all subgap current structures disappear and, more- 
over, the excess current at large voltage becomes neg- 
ative (deficiency current). The magnitude of the defi- 
ciency current is considerably large at <f> — tt. This effect 
results from the suppression of the Andreev reflections 
from the SiN and S2N interfaces due to the interference 
in the arms of thej-beam splitter (cf. a similar effect in 
NS interferometer^rQ). 

We conclude our discussion with some brief comments 
on the connection of the studied model to more realistic 
experimental devices based on 2D electron gas structures. 
In such devices, the normal region has a typical size which 
exceeds the superconducting coherence length and it con- 
tains several conducting electronic modes. Furthermore, 
the NS interfaces are not perfectly transparent. However, 
this effect is less important since the interface reflection 
coefficient in practice can be rather small (< 0.2). In- 
creasing the length of the arms of the interferometer will 
result in the appearance of new phase dependent cur- 
rent peaks due to increasing number of Andreev states 
in the SNS junction of the interferometer (cf. Ref. p^ . 
Similarly, increasing the length of the injection lead will 
result in the appearance of current resonances which do 
not depend on the phase difference, due to ^-independent 
superconducting bound states in this lead.Ej In inter- 
ferometers with multimode wires, the resonant current 
peaks will be smeared; however, one should still expect 
some phase dependent resonant feaiures similar to the 
case of diffusive NS interferometersQ. The interference 
effect leading to the large deficient current at large ap- 
plied voltage will survive in multimode junctions. 



VI. ACKNOWLEDGEMENTS 

The authors gratefully acknowledge support from NFR 
and KVA (Sweden) and NEDO (Japan). 



^ Th. Martin and R. Landauer, Phys. Rev. B 45, 1742 

(1992) . 

^ M. Biittiker, Phys. Rev. B 46, 12485 (1992). 

^ R.C. Liu, B. Odom, Y. Yamamoto, and S. Tarucha, Nature 

391, 263 (1998). 
* J. Torres and Th. Martin, Europhys. J. B 12, 319 (1999). 

H. Nakano and H. Takayanagi, Phys. Rev. B 47, 7986 

(1993) . 

^' F. W. J. Hekking and Yu. V. Nazarov, Phys. Rev. L 71, 
1625 (1993). 

A.F. Andreev, Sov. Phys. JETP 19, 1228 (1964). 
^ V.T. Petrashov et al., Phys. Rev. Lett. 70, 347 (1993). 
^ P.G.N, de Vegvar, et al., Phys. Rev. Lett. 73, 1416 (1994). 



11 



C.J. Lambert and R. Raimondi, J. Phys. Cond. Mat. 10, 

901, 1998 

" T.M. Klapwijk, G.E. Blonder and M. Tinkham, Physica 

B&C 109-110, 1657 (1982) 
12 J.M. Rowell and W.L. Feldmann, Phys. Rev. 172, 393 

(1968). 

" E.N. Bratus', V.S. Shumeiko and G. Wcndin, Phys. Rev. 

Lett. 74, 2110 (1995). 
^'^ A. Ingerman, G. Johansson, V. S. Shumeiko, and G. 

Wendin, Phys. Rev. B 64, 144504 (2001). 
1^ A phasc-dcpcndcnt singlo-particlc tunnel current in 4- 

terminal SINS interferometers with diffusive normal region 

was recently studied in: F.K. Wilhelm and A. A. Golubov, 

Phys. Rev. B 62, 5353 (2000). 
1^ J. Kutchinsky ct al., Phys. Rev. B 56, R2932 (1997). 
1^ H. Takayanagi, E. Toyoda, and T. Akazaki, Superlatt. and 

Microstr. 25, 993 (1999). 

A. Kadigrobov et al., Phys. Rev. B 52, R8662 (1995). 

i'"* P. Samuclsson, V.S. Shumeiko and G. Wcndin, Phys. Rev. 
B 56, R5763 (1997): Phys. Rev. B 62, 1319 (2000). 

20 L.F. Chang and P.F. Bagwell, Phys. Rev. B 55, 12678 
(1997). 

21 D. Averin and A. Bardas, Phys. Rev. Lett. 75, 1831 (1995). 

22 V.S. Shumeiko, E.N. Bratus', and G. Wendin, Low Temp. 
Phys. 23, 181 (1997). 

2^ G. Johansson, E.N. Bratus', V. S. Shumeiko, and G. 

Wendin, Superlatt. and Microstr. 25, 905 (1999). 
2* G. Johansson, E. Bratus', V. Shumeiko, and G. Wendin, 

Phys. Rev. B 60, 1382 (1999). 
2* P. G. de Gonnos and D. Saint-James, Phys. Lett. 4, 151 

(1963). 

2^ P.G de Gennes, Superconductivity in Metals and Alloys, 

Addison- Wesley Publishing Co., inc., 1989. 
2^^ J.R. Schrieffer and J.W. Wilkins, Phys. Rev. Lett. 10, 17 

(1963). 



12 



